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SPECTRAL APPROACH TO LINEAR PROGRAMMING BOUNDS ON CODES 


ALEXANDER BARG* * AND DMITRY NOGIN^ 


Abstract. We give new proofs of asymptotic upper bounds of coding theory obtained within the frame 
of Delsarte’s linear programming method. The proofs rely on the analysis of eigenvectors of some finite¬ 
dimensional operators related to orthogonal polynomials. The examples of the method considered in the paper 
include binary codes, binary constant-weight codes, spherical codes, and codes in the projective spaces. 


1. Introduction. Let A be a compact metric space with distance function d. A code C is a finite subset 
of X. Define fhe minimum disfance of C as d{C) = mmx^y£C,x^y y)- A variefy of mefric spaces fhaf 
arise from differenf applicafions include fhe binary Hamming space, fhe binary Johnson space, fhe sphere 
in M”, real and complex projective spaces, Grassmann manifolds, efc. Esfimafing fhe maximum size of fhe 
code wifh a given value of d is one of fhe main problems of coding fheory. Let M be the cardinality of 
C. A powerful technique to bound M above as a function of d{C) that is applicable in a wide class of 
metric spaces including all of the aforementioned examples is Delsarte’s linear programming method ^2!|. 
The first such examples to be considered were the binary Hamming space Hn = {0,1}"^ and the Johnson 
space J”’"' C Hn which is formed by all the vectors of Hn of Hamming weight w, with the distance given 
by the Hamming metric. The best currently known asymptotic estimates of the size of binary codes and 
binary constant weight codes were obtained in McEliece, Rodemich, Rumsey, Welch lUTl and are called the 
MRRW bounds. Shortly thereafter, Kabatiansky and Eevenshtein fT\ established an analogous bound for 
codes on the unit sphere in with Euclidean mehic and some related spaces. This paper also inhoduced a 
general approach to bounding the code size in distance-transitive metric spaces based on harmonic analysis 
of their isometry group. This approach was furthered in papers (U which also explored the limits of 
Delsarte’s method. 

In this paper we suggest a new proof method for linear programming upper bounds of coding theory. 
Our approach, which relies on the analysis of eigenvectors of some finite-dimensional operators related to 
orthogonal polynomials arguably makes some steps of the proofs conceptually more transparent then those 
previously known. We also consider some of the main examples mentioned above. The linear-algebraic 
ideas that we follow were introduced in a recent paper by Bachoc |T| in which a similar approach has been 
taken to establish an asymptotic bound for codes in the real Grassmann manifold. 

2. A bound on the code size. We assume that A is a distance-transitive space which means that its isometry 
group G acts doubly hansitively on ordered pairs of points at a given distance. In this case the zonal spher¬ 
ical kernels Ki (x, y) associated with irreducible regular representations of G depend only on the distance 
between x and y. In all the examples mentioned above, except for the Grassmann manifold, Aj(x,y) can 
be expressed as a univariate polynomial Pi{x) of degree i, where x = r((i) is some function of the distance 
d(x,y). 

Eet D be the (finite or infinite) set of the possible values of the distance in A. We will assume that 
r((i(x, y)) is a monotone function that sends ZJ to a segment [a, b], Eor instance, for the Hamming space, 
D = {0,1,... , n} and r can be taken the identity function. Eor the sphere D = [0, 2]. In this 
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case it is convenient to take T{d) = 1 — cP/2 to be the scalar product (x, y) = The invariant 

measure on G induces a measure d/r on [a, b]. For instance, for X = {0,1}”, the measure corresponds to 
the binomial probability distribution on {0,1,..., n}, so dfi = 1. We will assume that the last condition 
holds in general and normalize /r when this is not the case. 

The kernels iFj(x, y), i = 0,1,..., are positive semidefinite which means that fTj(x, y) > 0 for 

any finite set C C X. This property together with the fact that y) can be expressed as a polynomial 
of one variable gives rise to the following set of inequalities 

(1) pi(T(d(x,y))) > 0, 1 = 0,1,... 

x.yec 

called the Delsarte inequalities in coding theory. 

The function r can be chosen in such a way that the polynomials pi,i = 0,1,..., are orthogonal on 
[a,b] with respect to the scalar product {f,g) = f fgdp. Below we denote by V the space L 2 {dp) of 
square-integrable functions on [a, 6]. 

We will assume that the polynomials pi are orthonormal, i.e., \\pi\\‘^ = {pi,Pi) = 1- Note that this implies 
that pq = 1. Another assumption used below is that the product piPj for alH, j > 0 expands into the basis 
{pi} with nonnegative coefficients, i.e., 

(2) PiPj = qljPk idij > 0). 

k 

This property is again implied by the fact that the zonal spherical kernels are positive semidefinite, see B- 
Since the polynomials {pi} are orthogonal, they satisfy a three-term recurrence ifT^ of the form 

(3) xpk = akPk+i + PkPk + IkPk-i (/c = 0,1,... ;p_i = 0). 

Let Pi = epi , where e > 0 is some constant. We will write this recurrence in the form 

(4) PiPk = akPk+i + bkPk + CkPk-i, 

which follows from 0 upon noticing that Pi is a linear function. By ©, the coefficients ak,bk, Ck are 
nonnegative. 

Let C C X be a code of size M and distance d. Denote by A(C') = {r(d(x, y)), x, y e C, x / y} the 
set of values that the function r takes on the distances between distinct code points. Let tq = r(0). 

The main theorem of the linear-programming method asserts the following. 

Theorem 1. El El Let C C X be a code of size M. Let F{t) = YIILq FiPi{x) be a polynomial that satisfies 

(i) Fo > 0, Fi > 0,i = 1,2,, m; 

(ii) F{x) < Ofor X e A(C'). 

Then M < F{to)/Fq. 

The proof is obvious because on the one hand, by assumption (ii) 

Y F{T{d{x,y))) < MF{to)-, 

x,yec 

on the other hand, because of O, assumption (i) and the fact that po = 1, 

Y >FoM^. 

x,yeC i X,y 

This theorem is equivalent to a duality theorem for a linear programming problem whose variables are 
the coefficients of the distance distribution of the code C and whose constraints are given by the Delsarte 
inequalities. For this reason, estimates obtained from this theorem are called the linear programming bounds. 
Our objective in this section is to present a new method of obtaining bounds on M based on this theorem. 
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We shall use a generic notation bi,ai) for a tridiagonal matrix of the form 


Ak 


bo clq 0 0 ... 0 

Cl bi oi 0 ... 0 

C2 62 02 ... 0 

. Ofc-l 

0 0 . Cfc bk 


The largest eigenvalue of a square symmetric matrix M will be denoted by A ma v(M). 

Throughout the paper we use bold letters to denote operators acting on V and regular letters to denote 
their matrices in the basis {pi}. Let Vk be the space of polynomials of degree < k considered as a subspace 
of V. Let Efc be the orthogonal projection from V to Vk- Consider the operator 


Sfc = Efc o Pi : 14 ^ Vk, 


i.e., multiplication by Pi followed by projection on I4. The argument that follows relies on the fact that this 
operator is self-adjoint (with respect to the bilinear form (•,•)). Indeed, both multiplication by a function and 
the orthogonal projection are self-adjoint operators. Therefore, the matrix Sk = Ak{ci, bi,ai) is symmetric. 
In other words, 

ai = {Pipi,pi+i) = {pi,PiPi+i) = Cj+i. 

Ap X p matrix 4 > 0 (i.e., a matrix with nonnegative entries) is called irreducible if for any partition of 
the set of indices {1, 2,... ,p} into two disjoint subsets I and J, |/| -|- | J| = p, the matrix is 

nonzero (in other words, a directed graph G with vertices {1, 2,... ,p} and edges {i,j) whenever Aij > 0 
is strongly connected). For instance, the matrix Sk is nonnegative and irreducible. 

In the next lemma we collect the properties of irreducible matrices used below. 


Lemma 1. Let A>Obeapxp irreducible symmetric matrix. 

(a) Its largest eigenvalue Amax(^) A positive and has multiplicity one. There exists a vector y > 0 such 
that Ay = Amax(^)y- 

(b) X 

max (4) < maxi<j<p Ylj 

(c) For any y / 0, Amax(^) > 

(d) IfO < B < Afar some matrix B, or if B is a principal minor of A, then |Amax(f?)| < Aniax(^)- 


Here claims (a),(b),(d) form a part of the Perron-Frobenius theory (see, e.g., |!5l|), and claim (c) is obvious 
and holds true for any symmetric matrix. 

The suggested method for deriving upper bounds is based on the following theorem. 


Theorem 2. Let C C X be an (M, d) code and let pk = CLkPh+iifo)/Pkifo). Then 


M < 


4pfcPfc(ro) 

Pi(ro)-A 

max (Sk) 


far all k such that Aniax('S'fc-i) > Pi{x)forallx G A((7). 


Proof: Let g = Yli=i9iPi ^ ^k- Fix some p > 0 (its value to be chosen later). Consider the operator 
Tfc : 14 ^ 14 defined by 


(^) k9 ^k9 P9kPk, 

and let 9k be its largest eigenvalue. Recall that Tk is the matrix of this operator in the basis {pi}. (Tk is the 
same as Sk except that {Tk)k+i,k+i = {^k)k^i ^+1 “ P-) “shift” the matrix Tk by a multiple of 

the identity matrix I to make all of its elements nonnegative. For instance, we may consider Tk + Pkl > 0. 
Therefore, by Lemma ^d) we have 

'^max (Sk-i + pi) < 9k + p < -^max {Sk + pi), 
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whence we get 

(6) Aniax(*S'fc—l) < ^ ■^max(*S'fc)- 

Moreover, the eigenvalue 9^ is of multiplicity one. Denote by / = (/o, /i,..., fk) G 14 the eigenvector 
that corresponds to it. By (EJi we have 

Plf — (^kf “h PfkPk “h /fcflfcPfc+l) 


SO 


/ 


PPk “h 0,kPk+l 

Pi-Ok 


fk- 


Consider the polynomial F = {ppk + akPk+i)f- By Lemmata), / can be chosen to have positive coor¬ 
dinates. Therefore by the coefficients of the expansion of F into the basis {pi} are nonnegative. Next, 
if Amax(5'fc-i) > Pi{x) for X G A(C'), then by ® we have F{x) < 0 for x G A(C'), i.e., F{x) satisfies 
condifion (ii) of Theorem^ Since mulfiplicafion by / is a self-adjoinf operator, we compufe 


Fo = {{pPk + akPk+i)f, 1) = {pPk + a-kPk+iJ) = Pfk>0 

and 

P(- \ _ (PPkixo) + akPk+i{ro) f , ^ {pPkiro) + akPk+i{To))^ , 

Piiro)-Ok Pi(ro)-A^ax(5fc) 

provided fhaf Amax(5'fc) < n. Thus, 

F{to) jpPkjTo) + akPk+i{'ro)f 
Fo p(Pi(to)- -^max (-5^)) ■ 

The value of p minimizing fhe leff-hand side is p = pk- The claimed esfimafe is obfained by using fhe 
polynomial F = {pkPk + o-kPk+i)! in Theorem[2 I 

Remark 1. Nofe fhaf by Lemma^d), fhe {Amax('S'fc)} form a monotone increasing sequence. Therefore, 
fhe lasf condifion of fhe fheorem holds for all k greater than some value ko ■ 

Next let us estimate the largest eigenvalue of Su- 


Lemma 2. Let Oj+i > Oj, > 6^, i = 0,1,.... Then for aZZ s = 1,..., /c + 1, 

— (2(s I'jCLk—s+i T s-i-i) — Aniax('5fc) ^ o-fc—1 T max(c!,/j_]^ + bk—ij hk')- 
s 

Proof : By Lemma^b) 

Amax('S'A:) ^ max(afc_2 T bk—i T Ofc—i, Ufc—i T Zi^), 

hence the upper bound. On the other hand, take y = where t denotes transposition. Then by 

part (c) of the same lemma, 

^ z S—1 S—1 

-^max i^k) — (2 ^ ^ dk—p T ^ ^ Ok—p 

^ ^ p=l p=0 

Since we assumed that the coefficients ai , bi are monotone increasing on i this implies the lower bound, i 

Remark 2. In effect. Lemma |2l provides an estimate of the extremal zero of Pk+i- Indeed, consider the 
operator o x : Vfc ^ 14 - h is self-adjoint, so its matrix in the basis {pi} is tridiagonal symmetric 

and is given by Xk = Aki'^i, Pi,ai), where the elements ai,Pi,'yi are the coefficients in the three-term 
recurrence ©. 

It is well known (e.g., jhj]) that the spectrum of Xk coincides with the set of zeros of Pk+i- [A proof goes 
as follows: let pfc+i(A) = 0. Consider the action of on the polynomial / = pk+i/{X — x) £ Vk : 

Xf - Xkf = Xf- Efc(x/) = Efc((A - x)f) = FkPk+i = 0. 
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Conversely, if / G 14, / ^ 0, and 0 = A/ — X^/ = Efc((A — x)f), this implies that (A — x)f is a constant 
multiple of Pk+i- Therefore, pk+i{x) is proportional ^ to det(x/fc+i — X^).] Then the largest zero of 
Pk+i can be found as = Amax(Xfc) = max||j^||=i(Xfc?/, y), or more concretely as 

^ k k-1 

xt+i = max V Pivf + 2 V amyi+i 

11^11=1 I 7^0 7^0 

This formula was first published in Gi p.580] with a different proof. 

We note that the relation between the extremal zero of Pk+i and the largest eigenvalue Amax(Xfc) makes 
the task of finding the zero computationally much easier that the direct approach because of the existence 
of very efficient iterative algorithms for the symmetric eigenvalue problem. This property is helpful for 
computing linear programming bounds on codes such as the bounds considered in the next section and other 
similar results for codes of moderate or even large length (on the order of several thousands). 


3. Examples. In this section we consider a few examples of interest to coding theory. 


3.1. Binary codes. Let X = {0,1}” be the binary Hamming space. It is known j2l|7l that the polynomials pi 
are given by the (normalized) Krawtchouk polynomials {Kk (x), A: = 0,1,..., n}. We have p{i) = , 

so the bilinear form can be written as (/, g) = m(*)/(*) 5'(0' Let C be a binary code of length n, size 

M and minimum Hamming distance d = d{C). We choose r(/c) = A: to obtain A(C') C {d, d + 1,..., re}. 
This inclusion may be proper depending on the code C, but we will ignore this and assume that A(C') = 
{d, d + 1,..., re} since this assumption can only relax the linear programming bound on M. 

The polynomials Kk satisfy a three-term recurrence relation O 

(7) 2xKk{x) = -a/( re - k){k + l)Kk+i{x) + nKk{x) - {n - k + l)kKk-i{x), 

Ko = 1, ki{x)kj{x) = Yjk dijKkix) with q^j > 0, and 


( 8 ) 



Choose in © Pi 
i/{i + l)(re - i),i = 

r 0 


Sk = 


0 


0 


= \kipi = re — 2x. From Q we then obtain Sk = Ak{ai-i,0,ai), where a* = 
0, 1 ,..., or more explicitly, 

y/n 0 ... ... 0 

0 y/2{n - 1) ... ... 0 

y/2{n - 1) 0 ... ... 0 

0 y/ {k — l){n — k + 2) 0 

... {k — l)(re — A; + 2) 0 y/ k{n — k + 1) 

0 ... 0 ■\/k{n — k + 1) 0 


The monotonicity assumption of Lemma|2lclearly holds because ak > Ok-i as long as A: < re/2. Therefore 
for the largest eigenvalue of Sk we obtain the following estimate: 


-—•a/(A; - s + 2)(re - A: + s - 1) < Amax(*S'fc) < 2y/k{n - k + 1). 

s 

Letting n ^ oo, s ^ oo, s = o{n), we obtain the exact asymptotic behavior of the main term: 
(9) hm = 2^t{1-t). 

n—>oo, kjn^T Tl 


1 


The coefficient equals aoai 


Q.k -1 and can be found recursively from and the equality po = 1. 
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Since tq = 0 and = n — k, the bound of Theorem|2ltakes the form 


( 10 ) 


M < 


4(n — k) /n 


^ -^max V 

for all k such that Amax(5'fc-i) > Pi{d) = n — 2d. This estimate together with ® leads to the following 
asymptotic result (the asymptotic MRRW bound for binary codes illl l: 


— logM < h {^/2 - \/5{l - (5))(1 + o(l)). 
n 


Here h{x) = x log 2 x — (l — x) log 2 (l — x) is the binary entropy function. Indeed, let lim - = 6 and assume 


that 6 < 1/2. We need to choose k so that > (1 — 25)(1 + o(l)) as n —> oo. In the limit, this amounts 
to taking r that satisfies 2^J r(l — r) > 1 — 25, or r > 1/2 — 5(1 — 5). The result now follows by the 

Stirling approximation. 

Remark 3. Specializing Remark 2 to the case at hand, we observe from 0 that 


Xk = V2(n4+i - Sk) = y 2 Aki-\/iin - i + l),n,+ l){n - i)). 
Therefore we obtain the following expression for the largest root of Kk+i ■ 


k-l 


^k+1 = ? + “F - i). 

^ V =1 


i=0 

This result is originally due to ^\. Although more accurate estimates of the extremal zeros are available in 
the literature fTOl l^. our Lemma|2l suffices to compute the correct value of the main term. 

Remark 4. The bound (IlOb is close to the previously known estimates obtained within the frame of 
Delsarte’s method. In particular, Levenshtein (8l|l0l constructed a sequence of polynomials that are optimal 
in the Delsarte problem (with some qualifiers). His results imply that the above estimate does not improve 
the known bounds on M. The result of ca is also of the form similar to (flfll . 

Remarks 2-4, modified appropriately, apply also to the other examples in this section. 

3.2. Constant-weight codes. Now let X C J”’"' the binary Johnson space, i.e., the set of vectors in 
{0,1}*^ of Hamming weight w. We take d to be the Hamming metric so that D = {0, 2,... , 2w} and 
put T{d) = d/2. The relevant family of orthogonal polynomials is given by the Hahn polynomials H^x) 

{w\(n—w\ 

— ^ according 


|2!|. They are orthogonal on t{D) = {0,1,... , u)} with respect to the weight Pj{i) = 
to f HkHmdpj = {y)dkm and satisfy a three-term recurrence 


0 


(11) {k + l){w — k){n — w — k){n — 2k + 2){n — 2k + 2>)Hk+i{x) = 

{n — 2k — l)(n — 2k + 3)[(n + 2)w{n — w) — nk{n — k + 1) — {n — 2k){n — 2k + 2)x]Hk{x) 

— {n — 2k — l){n — 2k){w — k + l){n — w — k + l){n — k + 2)Hk-i{x). 

Note that Y17=i = 1- Let us normalize by setting (fc)) As above, we have 

W 

Hyx)Hjix) = dlAix) (qlj > 0) 


k=0 


and 


Let us take 


^A:(0) = 


in — 2k + 1 


n — k + 1 
Pi{x) = {n- l)-^/^Hi{x) = 1 - 


nx 


w{n — w) 
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Let us write out the matrix of the operator = E^. o Pi in the orthonormal basis. We have = 
Ak{ai-i, bi,ai), where the matrix elements can be computed from (fTTl . We obtain 

{i + l)(n - z + 1) 

(n-2z + l)(n-2z- 1)’ 

, (n - 2u;)^i(n - i + 1) 

bi = -^-77-77-7, z > 0. 

w{n — w){n — 2z)(n — 2z + 2) 

Let C C J”’"' be a code of size M and distance 2d. Let us apply Theorem |2to bounding M as a function 
of d. We have tq = 0,Hq = 1, 

Hk+i{0) n{w — k){n — w — k){n — k + 1) 

Pk CLk w{n - w){n - 2k){n - 2k + 1) 

and A(C) = {0,1,... , d}. Thus, we obtain the following estimate. 


di — 


n{w — i){n — w — i) 
w{n — w){n — 2i) 


Theorem 3. 

‘in{w — k){n — w — k) / n\ 

“ (1 - ^max(Sk))w(n - w)(n - 2k) [kj 
for all k such that Amax(*S'fc-i) > 1 - ■ 

Let us find the minimum k that satisfies fhe required condition. Firsf we use Lemma |2 to compute the 
asymptotic behavior of Amax(5'fc). 


Lemma 3. 


lim 

n—>co 

ifi/n—, A:/?!—»■ r 


'^max (SO 


uj{l - u;)(l + 2y^r(l - r)) 


Proof: Note that for the upper bound in Lemma|2lis suffices fo prove fhaf fhe value ai + bi + aj_i grows 
on z. Letting a = ^, we compufe 


o-i-l + bi p Oi 


2{lo — a)(l — io — a)^a{l — a) + (1 — 2w)^a(l 
u;(l — a;)(l — 2a)^ 


-(1 + 0 ( 1 )). 


^ 2u;(l-u;)^Q((l-a) + 0((l-a) ^ 

a;(l — + 2i/a{l — a)) 

The main ferm on fhe righf-hand side of fhe lasf expression is a growing function of a. Indeed, a{l — a) 
grows on a for a < 1 / 2 , so we only need fo check fhaf fhe function f(2u;(l — u)) + t)/{l + 2t) increases on 
ffor 0 P t P ^^2 which IS sfraighfforward. Thus we put z — k — 1 and obtain for Aj^ax(*^A:) upper bound 
of the form claimed. Lemma|2lalso implies a matching lower bound. Namely, from its proof, 

Amax('5fc) ^ — ^2 'y ^ CLk—p p ^ ^ bk—p^ {s — 1, . . . , A) T 1). 

p=l p=0 

For large values of the parameters, we can write 

Amax(>S'fc) ^ T bk—s+1 T Ofc—s+l)(l T o(l))' 

The proof is completed by letting s ^ oo, s = o{n). I 

Let us use this lemma in Theorem |3l Assume that n ^ oo,d = 6n. The condition on k in this theorem 
will be fulfilled for any k = r/n fhaf safisfies 

> 1 - 

uj{l — ixi)(l + 2y^r(l — r)) uj{l — uj) 
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or 

^ ^ {uj - t){1 - U) - t) 

1 + 2^t{1 - t) 

We conclude that Theorem|3limplies the following estimate for an (n, M, 2dn) code C C J”’"' (the asymp¬ 
totic MRRW bound for constant weight codes fldll '): 

- log M < /i(r)(l -h o(l)), 
n 

where 5 = {u) — t){ 1 — uj — t) /{1 + 2y^r(l — r)). 


3.3. Spherical codes. Consider codes on the unit sphere 5” ^ in W^. The polynomials pi in this case belong 
to the family of Gegenbauer polynomials Cfc(x) fT^ pp.SOff]. We have 

1 ^_3 M-t-i-S-j 

Ci(x)Cj(x)(l — X^)^~ dx = —- - - LOndi j, 

^ j 3 \ ) n + 2i-2 

where ujn = n„,- 2 r 2 ^n- ^; particular for i = j = 0, /i^(l — dx = lOnjin — 2). We also have 

ct(i)=ry’). 

Normalizing the measure, we obtain dp{x) = ^^(1 — The normalized Gegenbauer poly¬ 

nomials are then given by 

- 2 

+k-3 
k 

The polynomials Ck satisfy a three-term recurrence of the form 

xCk{x) = ttkCk+iix) + ak-iCk-iix), 

where a* = i = 0,..., and C-i = 0, Cq = 1. Further, CiCj = where q’^j > 0 

and 

n + 2k — 2 /n + k — 3\ 
n — 2 \ k )' 

Let C{n, M, t) denote a code in which the angle between any two distinct vectors Xj, xj satisfies cos(xj, x^ ) 
< t. As remarked above, we take r(d) = 1 — d‘^/2. We have D = [0,2],r(-D) = [—1,1],A(C') C 
[—1, f], To = 1. Choose Pi{x) = = x, then the matrix Sk has the form Ak{ai-i,0, Oj), so 

Ca:+i(1) n + k-2 

^ Cfc(l) n + 2A:-2 

From Theorem |2l we obtain 




^ _ I n + 2k 


Theorem 4. 

( 12 ) 


M < 


1 ^max('5fc) 

for all k such that Amax(>S'fc-i) > t. 

This coincides with the original bound of fj]. 

Lemma 4. For any s = 2,..., k 

2{s — 1) I {n + k — s — l){k — s + 2) 
s y {n + 2k — 2s + 2){n + 2k — 2s) 


n + k — 2 
k 


< Amax(5fc) < 2. 


{n + k — 3)k 


{n + 2k — 2)(n -|- 2A; — 4)' 
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In particular, 


lim 


'^max jSk) ^ ,^y/p(l+p) 


n 


n—^oo, —— 

’ n ^ 

Proof: We only need to check that a, > Oj+i. For n > 5, 


l + 2p 


a, - a,_i = 


(n — 2)(n — 4) 


> 0 , 


(n + 2i)(n + 2i-2)(n + 2z-4) 
so Oj is an increasing function of i. The inequalities in the claim now follow directly from Lemma|2 Letting 
s oo, s = o(n) and taking the limit gives the asymptotic behavior of Amax( 5 'fc). I 

Theorem|4]and Lemma |4] together enable us to recover the asymptotic bound of |7i. Namely, using the 
Stirling approximation we obtain 

1 


n 


logM < ((1 + p)log(l + p) - plogp)(l + o(l)) 


under the condition t < Amax(-^fe-i) which in the limit of n ^ oo, - —> p translates into p > 


i-vT^ 

2vT^ ■ 


3.4. Codes in projective spaces. A class of spaces related to the real sphere is given by the projective spaces 
where L = M or C of H. The zonal spherical functions in these spaces are given by the Jacobi 
polynomials ffTl . where a = a{n — 1) — 1, /J = u — 1, and a = 1 / 2 ,1, 2, respectively. 

The polynomials satisfy 

f P°‘’^(x)P°’^(x)(l - x)“(l + x)^dx = _ 2^+^+ (fc + a)!(fc + /J)! _ 

X) [L + x) dx ^2k + a + (3 + l)kl{k + a + P)r^' 

, , fk + a 
Pk{l) = 

\ a 

where by definition x! = r(x + 1). The coefficients of three-term recurrence Q have the form 

_ 2{k + 1) (k + q; + /J + 1) _ pf — of 

“ {2k + a +p + l){2k + a +p + 2y ~ {2k + a + P){2k + a + P + 2)' 

2{k + a){k + P) 

{2k + a + P){2k + a + P + 1) 

Define fhe bilinear form on V by (/, g) = fgdp, where 

Then the squared norm of Pk is equal to 


\P 


,«,/3||2 _ (q + /? + 1)(« + Pf{k + a.)\{k + P)\ 
{2k + a + P + l)a\p\k\{k + a + P)\ 


Denote by Pk = P^'^ j 
We will take in @ 


the normalized Jacobi polynomials. 


1 


= = 

then the coefficients of the recurrence are found to be 

a + /3 + 2 


Pi{x) = P^ ’’^{x) = -{{a + P + 2)x + a- P), 


a-k = 


2k + a + P + 2\ 

h = 


(/c + q; + f){k + /3 + l)(/u + l){k + cr + /3 + 1) 
{2k + a + P + 3){2k + a + P + 1) 

2{a — P)k{k + a + P + 1) 


{2k + a + P){2k + a + /3 + 2)' 
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and Ck = ak-i. 

Let C C X be a code of size M in which | (xj, Xj) | < t for any two distinct vectors x*, xj. We have 
D = [0, -v/2], so choosing r(d) = 2(1 — (fi j2)‘^ — 1 we obtain t{D) = [—1,1], A(C') C [—1, 2t^ — 1], We 
compute 

p2,,, ^ (2t + a+ /?+!) (TiCT'’) 

’ Q + /3+1 (W) 

(a + /3 + 2) (A) + a + 1)(A: + o; + /? + 1) 

^^(1) {‘2k + a + /? + 1)(2A; + a + /3 + 2) 

Using these expressions in Theorem|2we obtain 


Theorem 5. 

^ 4(a + /3 + 2){k + g + 1) 

{2k + a + /? + 2)(1 — Aniax('S'fc) 

Let us use LemmaEJto derive the asymptotic behavior of Amax(*S'fc) as A: —> oo, a = ak, j3 = bk, a > 
0, 6 > 0. We obtain 


Amax(5'fc) ^ 2((a + b )^(a + 1)(6 + l)(a + 6 + 1) + {a — b){a + b + 1)) 
k (a + 6 + 2)^ 

The condition for Theorem |2l to be applicable is 
(13) Xra^ASk) > Pi{2t^ -I) = {a + P + 2)^ -(3-1. 

Forinstance, let us derive a bound for the case X = Letting A; = snl2,a = {n — 3)/2, (3 = — 

we obtain a = Vs, 6 = 0, 

Amax(<S'fc) 4(1 + a) 
k ' (1 + 2s)2' 

Therefore, for large values of the parameters condition (d becomes 

4(1 + s) _ V 

(1 + 2s)2 “ 7’ 


or s = V2((l/Vl — P) — !)• From Theorem|2lwe obtain the asymptotic bound of 0 on the code size: 

— log M < (1 + s) log(l + s) — s log s. 
n 

In a similar way we can recover the asymptotic bounds of 0 in the other cases mentioned. 

The method presented is a linear-algebraic alternative to the analytic methods of El III 113 It is equiva¬ 
lent to them in the sense that it gives the same asymptotic results, although for finite parameters the bounds 
derived by these two approaches generally do not coincide. 
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